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Abstract 

We consider scenarios in which the lightest sparticle (LSP) is the gravitino and the next-to- 
lightest sparticle (NLSP) is a metastable stau. We examine the production of stau pairs in e + e~ 
annihilation at ILC and CLIC energies. In addition to three minimal supergravity (mSUGRA) 
benchmark scenarios proposed previously, we consider a new high-mass scenario in which effects 
catalyzed by stau bound states yield abundances of 6 ' 7 Li that fit the astrophysical data better 
than standard Big-Bang nucleosynthesis. This scenario could be probed only at CLIC energies. 
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In each scenario, we show how the stau mixing angle may be determined from measurements of 



Qh the total stau pair-production cross sections with polarized beams, and of the tau polarization 



in stau decays. Using realistic ILC and CLIC luminosity spectra, we find for each scenario the 
center-of-mass energy that maximizes the number of staus with (3"f < 0.4, that may be trapped 
in a generic detector. The dominant sources of such slow-moving staus are generically the pair 
production and cascade decays of heavier sparticles with higher thresholds, and the optimal 
center-of-mass energy is typically considerably beyond 2mf 1 ■ 
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I. INTRODUCTION 



Supersymmetry (SUSY) is one of the most interesting possibilities offered by quantum 
field theory whose relevance to nature has not yet been established PJ. In addition to 
relating fermions to bosons, stabilizing the mass scale of the weak interactions and (if Im- 
parity is conserved) providing a plausible candidate for cosmological dark matter (DM), 
supersymmetry also provides a framework for the unification of particle physics and gravity. 
The latter is governed by the Planck scale, M P ~ 10 19 GeV, which is close to the energy 
where the gravitational interactions may become comparable in magnitude to the gauge 
interactions. 

If supersymmetry is to be combined with gravity, it must be a local symmetry, and the 
resulting theory is called supergravity (SUGRA) In models of spontaneously- 

broken supergravity, the spin-3/2 superpartner of the graviton, namely the gravitino G, 
acquires a mass by absorbing the goldstino fermion of spontaneously-broken supersym- 
metry via a super-Higgs mechanism {fj,Lzl]. The gravitino is a very plausible candidate for the 
DM, since the presence of a massive gravitino is an inevitable prediction of SUGRA. More- 
over, in models of gauge-mediated SUSY breaking [8| the gravitino is necessarily the LSP. 
Furthermore, there are also generic regions of the parameter space of minimal supergravity 
(mSUGRA) models where the gravitino is the LSP. 

In mSUGRA the soft supersymmetry-breaking parameters take a particularly simple form 
at the Planck scale [5j, in which the soft supersymmetry-breaking scalar masses-squared rn|, 
gaugino masses my 2 and trilinear parameters A are each flavor diagonal and universal [9||. 
Renormalization-group evolution (RGE) is then used to derive the effective values of the 
supersymmetric parameters at the low-energy (electroweak) scale. Five parameters fixed at 
the gauge coupling unification scale, m , the gaugino mass m^, A , tan/3 and the sign of 
(/i) are then related to the mass parameters at the scale of electroweak symmetry breaking 
by the RGE. Moreover, mSUGRA requires that the gravitino mass mg = mo at the Planck 
scale, and imposes a relation between trilinear and bilinear soft supersymmetry-breaking 
parameters that may be used to fix tan (3 in terms of the other mSUGRA parameters. The 
mSUGRA parameter space offers both neutralino dark matter and gravitino dark matter as 
generic possibilities. 

In general, it would be difficult to detect gravitino DM directly, since the gravitino's cou- 
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plings are suppressed by inverse powers of the Planck scale. However, evidence for gravitino 
DM may be obtained indirectly in collider experiments, since in such a case the next-to- 
lightest sparticle (NLSP) has a long lifetime, as it decays into the gravitino through very 
weak interactions suppressed by the Planck scale. One of the most plausible candidates for 
the NLSP is the lighter scalar superpartner of the tau lepton, called the stau, T\ . If universal 
supersymmetry-breaking soft supersymmetry-breaking masses for squarks and sleptons are 
assumed, as in mSUGRA, renormalization effects and mixing between the spartners of the 
left- and right-handed sleptons cause the t\ to become the NLSP. 

If R-parity is conserved, sparticles must be pair-produced at colliders with sufficiently 
high energies, and heavier sparticles will decay into other sparticles, until the decay cas- 
cade terminates with the NLSP, which is metastable in gravitino DM scenarios based on 
mSUGRA and decays only long after being produced. Collider experiments offer possibilities 
to detect and measure metastable NLSPs such as the staus in mSUGRA, and to use their 
decays to study indirectly the properties of gravitinos, which cannot be observed directly in 
astrophysical experiments. 

There are several experimental constraints on the mSUGRA parameter space: the LSP 
must be electrically and colour-neutral, and one must satisfy the limit from LEP 2 on the 
stau mass: > 85.2 GeV obtained from the pair production of staus, as well as limits on 
other sparticles, the mass of the lightest Higgs scalar, the DM density, the branching ratio for 
radiative B decay and the muon anomalous magnetic moment. An up-to-date compilation of 
these experimental and phenomenological constraints on the mSUGRA model can be found 
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Astrophysics and cosmology also constrain the properties of metastable particles such 



as the staus in gravitino DM scenarios 



12| . in particular via the comparison between the 



calculated and observed abundances of light elements, assuming the cosmological baryon 
density derived from observations of the cosmic microwave background (CMB) . In addition 
to the hadronic and electromagnetic effects due to secondary particles produced in the 
showers induced by stau decay, one must also consider the catalysis effects of stau bound 



states 



13| . In certain regions of mSUGRA parameter space, these may even bring the 



calculations of Big-Bang nucleosynthesis (BBN) into better agreement with the observed 
abundances of 6 ' 7 Li [141. A number of benchmark scenarios in which the stau decay showers 

n 

do not destroy the successful results of BBN were proposed in [16fl. These did not take 
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bound-state effects into account, but these may be included in related models, as we discuss 
later. 

The phenomenology of stau NLSPs in both proton-proton and e + e~ colliders has been 
considered in a number of previous papers. In particular, it has been shown that the stau 
mass could be measured quite accurately by the ATLAS and CMS experiments at the 
LHC, using their tracking systems and time-of- flight (ToF) information from Resistive Plate 
Chamber (RPCs). Previous studies show that the polarization of r lepton could be measured 
via the energy distribution of the decay products of the polarized r lepton at future lepton 



colliders 



17 1 . There have also been some phenomenological investigations of long-lived staus 



at the ILC 



cQ, 
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In this work, we study the pair production of staus in e + e~ annihilation at the future 
linear colliders ILC {22] and CLIC 23], assuming that the LSP is the gravitino and the 
NLSP is the lighter stau. In addition to the three benchmark scenarios proposed earlier, 
we also propose a new benchmark scenario chosen to improve the agreement of Big-Bang 
nucleosynthesis calculations with the observed 6,7 Li abundances after bound-state catalysis 
effects are included. This point has a relatively heavy stau that could only be observed at 
CLIC. Following a general discussion of the signatures and detection of staus in the linear- 
collider environment, we explore the possible measurements of model parameters such as the 
stau mixing angle that could be made at the ILC or CLIC. We show that beam polarization 
and/or a measurement of the final-state r polarization would be very useful in this respect. 
We also discuss the choice of the center-of-mass energy that would maximize the cross section 
for producing slow-moving staus with /J7 < 0.4, that would be trapped in a typical linear- 
collider detector. The subsequent stau decays would be optimal for some measurements 
of gravitino parameters. We show that most such slow-moving staus are produced in the 
cascade decays of heavier sparticles, so that the optimal center-of-mass energies for trapping 
staus are typically considerably higher than the stau pair-production threshold. 



II. THE MODEL PARAMETERS 



In the mSUGRA framework used here, \i is the supersymmetric Higgs(ino) mass param- 
eter, we denote by m the universal soft SUSY-breaking contribution to the masses of all 
scalars (including the Higgs fields), is the universal SUSY-breaking gaugino mass, and 
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Aq is the universal SUSY-breaking factor in the soft trilinear scalar interactions. In addition, 
there is a bilinear soft SUSY-breaking parameter B , which appears in the effective Higgs 
potential, and is related by mSUGRA to the trilinear parameter: B = A — m . All these 
soft SUSY-breaking parameters are defined at the scale of grand unification. In general 
studies of the MSSM, tan/?, the ratio of the vacuum expectation values (vevs) of the two 
Higgs doublets at the weak scale, is taken as a free parameter. However, in mSUGRA it 
is determined dynamically once the relation between Aq and Bq is imposed. This param- 
eter set also determines the gravitino mass, since mg = m in mSUGRA, as well as the 
other sparticle masses. For calculations of sparticle mass spectra and other properties in 
the models we study here, we use the updated ISASUGRA programme, which is part of the 



ISAJET package 



24j. This programme is convenient for event generation with PYTHIA 



for various colliders. 

A selection of benchmark mSUGRA points consistent with present data from particle 



physics and BBN constraints on the decays of metastable NLSPs was proposed in [14j. We 
use a subset of these benchmark scenarios to study mSUGRA phenomenology at e + e~ col- 
liders, namely the points denoted (e, C,T]), in which the NLSP is the lighter stau, f±. The 
fi lifetime in these models ranges from ~ 10 4 s to ~ 3 x 10 6 s, so fast-moving charged NL- 
SPs would be indistinguishable from massive stable particles, in a generic collider detector. 
However, it might be possible to observe the decays into gravitinos of any slow-moving tlS 
that lose sufficient energy to stop in the detector. 

The mSUGRA benchmark scenarios were formulated assuming the representative value 
Aq = (3 — V3)m found in the simplest Polonyi model of supersymmetry breaking in 
mSUGRA 0| . The allowed region of the gravitino DM parameter space may be displayed 
in the (mi/ 2 , m ) plane, as shown in the left panel of Fig. [Qfor the Polonyi value of A . The 
theoretical and phenomenological constraints are displayed, together with the astrophysical 
constraints on NLSP decays. The combined effect of these constraints is to allow a wedge in 
the (mi/ 2 ,m ) plane at relatively low values of mo. The benchmark point e was chosen near 
the apex of this wedge, whereas the points r\ and ( were chosen at larger values of mi/ 2 , 
and hence more challenging for the LHC and other colliders [if]]. The mSUGRA parameters 
specifying these models and some of the corresponding sparticle masses are shown in Table HI 
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Figure 1: The left (first) panel shows the (mo, m^) plane for mSUGRA with a gravitino LSP and 
stau NLSP, assuming Aq = 3 — y/3, with the regions allowed by astrophysical constraints on stau 
decays shaded light (yellow). The values of tan (3 are fixed by the vacuum conditions with the 
values indicated along the solid black contours, and the three benchmark points e, r\ and £ are shown 
by red + signs. The right (second) panel shows the (mo,mu2) plane in the CMSSM, assuming fixed 
tan/3 = 10, with the regions allowed in the presence of cosmological bound-state effects and stau 
decays ifm^ = 0.2mo shaded light (yellow). In the darker (pink) shaded region, bound-state effects 
yield 6,7 Li abundances in the ranges favoured by astrophysical observations. The new benchmark 



point 9 is located in this 6,7 Li-friendly region \l. 



as calculated using ISA JET 7.74 1 



Recently, bound-state effects on light-element abundances in models with metastable 

nn 

stau NLSPs have been studied in several variants of the MSSM [1J, ll5|, including mSUGRA 
and some other models in which its relations between and m and between A and B 
are relaxed, but where universality is maintained for the soft SUSY-breaking parameters, 
a framework we term the constrained MSSM (CMSSM). The bound-state effects restrict 
greatly the allowed regions of parameter space, but one example was found within the 
CMSSM of an (mi/ 2 ,m ) plane with ttiq = 0.2m and tan/3 = 10, that has a wedge compat- 



[1] In addition, we note that the en — fx mass difference for point 9 is just 0(100 MeV). 
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Table I: The four benchmark points with gravitino DM and a f\ NLSP that are considered in this 
paper. The quoted slepton, squark, gluino and the neutralino mass spectra (in GeV) have been 
calculated with ISA JET 7.74, using mt = 175 GeV. 
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ible with BBN and also a small region where the BBN predictions for the ' Li abundances 
could be brought into better agreement with astrophysical observations. We propose and 
study here a new benchmark point in this 6 ' 7 Li-friendly region, denoted by 9, with the pa- 
rameter choices shown in the last column of Table [B Also shown in this column are some 
sparticle masses: like the rest of the 6,7 Li-friendly region, benchmark 9 features relatively 
heavy sparticles beyond the reach of either the LHC or the ILC, but within the kinematic 
reach of CLIC. 

In this connection we note that for benchmark 9 the direct stau pair production cross 
section is 2.7 x KT 6 pb at the LHC with y/s = U TeV and 2.8 x 1(T 5 pb at the DLHC 
with s/s = 28 TeV. The plausible integrated luminosity of 10 5 pb -1 for the LHC would be 
insufficient to reach this point, and it would be very marginal at the SLHC with 10 6 pb -1 . 
The only alternative to CLIC would be the DLHC with a very high integrated luminosity. 

The bound-state effects depend sensitively on the stau lifetime, being much more im- 
portant for longer lifetimes. This is why they would be unacceptable for the the stau in 
benchmark scenarios e, rj and £, if it had the lifetime predicted by mSUGRA. However, 
as discussed in the following Section, the stau lifetime depends sensitively on the mass as- 
sumed for the gravitino. Hence, a CMSSM variant with similar masses for the masses of 
the spartners of the Standard Model particles but a smaller gravitino mass would not be 
excluded a priori on the basis of bound-state effects. For this reason, and since the previous 
benchmarks e, r] and ( represent to a great extent the range of possibilities open to the 
ILC (limits on metastable particles exclude a stau much lighter than at benchmark e, and 
sparticles would be inaccessible to the LHC or ILC if they weighed much more than at points 
rj and (), we continue to include all these points in the subsequent discussion, though our 
primary interest will be in the new 6,7 Li-friendly point 9. 

III. THE n NLSP DECAY AND LIFETIME 

Since we assume the conservation of R-parity, the gravitino must be produced at the 
end of any decay chain initiated by the decay of a heavy unstable supersymmetric particle. 
As already remarked, in models with universal scalar soft supersymmetry-breaking masses, 
the lightest slepton is generically the lighter stau mass eigenstate, the r 1? as a result of 
renormalization-group effects and left-right sfermion mixing. In large regions of mSUGRA 
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parameter space, the T\ is the NLSP, and is the penultimate sparticle in (essentially) every 
decay chain. 

The T\ is in general a linear combination of tl and Tr, which are the superpartners of the 
left- and right-handed r leptons Tl and Tr, respectively. In general the stau mass eigenstates 
are 

ti = tl cos 6f + T R sin9f, (1) 
t 2 = —tl sin Of + t r cos Of, (2) 

where Of is the stau mixing angle, which is given by 

-m T (A T - fi tan (3) 
cos Of = , (3) 

J {m% L - m\) + m T (A T - /xtan/?) 2 

where | cos f \ > 1 / \/2 if m? L < rrif R and | cos 6?\ < 1/ V2 if rrif R < m^ L . 

The interactions of the stau states tl,r with the gravitino G and tau lepton r are described 



by the Lagrangian 



26] 



Lrrd V2M F 



(p v T L )*(fi'7 lt P L T + {D u tl)tPl1»1 v G» 



+ {D v t r )*G j u ^P r t + {D v t r )tP r1 ^G» , (4) 



where D u = (d u + ig e A u ), with A v denoting the gauge boson and Mp = (8tvGn) l ' 2 = 
2.436 x 10 18 GeV is the reduced Planck mass, with the Newton constant G N = 6.707 x 1(T 39 
GeV" 2 . 

The stau decay rate is dominated by the two-body decay into tau and gravitino (ri — > tG), 
and the decay width of the T\ is given by 

_ [m\ - m| - m 2 T ) 4 ( 4m|m^ \ 3/2 

T1 ~ 487rM2m|m| ^ (ml -m\- m\Y J ' { } 

where m^, and m T are the masses of the stau Ti, gravitino G and tau lepton r, respec- 
tively. As illustrative examples, if we take = 150 GeV (similar to the mass at benchmark 
point e), we find a lifetime of T- 1 = 7762 s (25.9 years) for a gravitino mass of mg = 1 
GeV (100 GeV), respectively. The lifetime and decay width of the stau as functions of the 
gravitino mass for various stau mass values are presented in Fig. El 

Neglecting the tau lepton mass, the mean stau decay length, obtained from the relation 
L = Pytf, is shown in Fig. [3l For instance, taking /3j — 1 and = 1 MeV would imply a 
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Figure 2: The stau NLSP lifetime and decay width as functions of the gravitino mass, for various 
choices of stau mass. 

mean decay length of L = 2.33 x 10 6 m (L = 109.6 m) for m fl = 150 GeV (m fl = 1100 GeV), 
respectively. Only for < 1 KeV might a typical stau with = 0(1) decay within a 
generic collider detector. In this case the light gravitino would constitute either warm or even 
hot dark matter, which is disfavoured by the modelling of cosmological structure formation. 

However, some fraction of the staus may be produced with velocities sufficiently low 
that strong ionization energy loss causes them to stop within the detector. In such a case, 
interesting information may be extracted from the stau decays. For example, the mass of the 
gravitino produced in the decay could in principle be inferred kinematically from the relation 
m~ = m| a + ml — 2m T E T , unless it is very small. Also, we note that the polarization of the 
tau produced in the stau decay must coincide with the combination of tl,r contained in the 
fi, offering in principle the possibility of determining the stau mixing angle 9f by examining 
the kinematical distributions in the subsequent tau decays, as well as in the polarization 
dependence of stau pair production. We discuss later the prospects for measuring 9? using 
the stau pair-production cross section in e + e~ collisions. 

The point e is near the tip of the wedge in Fig. [Q and the stau NLSP fi has a mass of 
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Figure 3: The mean decay length of the stau NLSP as a function of its mass, for various choices of 
the gravitino mass. 

m f = 154 GeV and a lifetime of t f = 2.9 x 10 6 s in mSUGRA 2 . The benchmark point ( 
is close to the upper edge of the wedge, with m = 100 GeV and tan/5 ~ 21.5, and the fi 
has a mass of 346 GeV and a lifetime of 1.7 x 10 6 s in mSUGRA. The benchmark point 7] 
is close to the lower edge of the wedge in Fig. CD, with m = 20 GeV and tan/3 = 23.7. Here 
the f\ has a mass of rrif = 327 GeV and a smaller lifetime of 6.4 x 10 4 s in mSUGRA. The 
new benchmark point 6 introduced here with a much larger value of m^, corresponding 
to rrif = 1140 GeV and tf = 1.4 x 10 3 s. The stau mixing angle Of takes the following 
values at these points: 77.2°, 82.4°, 81.8° and 89.0° for the points e, (, r\ and 9, respectively. 
Thus, each of these scenarios predicts that the f\ should be mainly right-handed. Then, the 
right-handed polarization of the electron beam enhances the f^f{ signal. We expect this 
feature to be quite general, and more marked for heavier staus. Correspondingly, the tau 
produced in stau decay should also be mainly right-handed. 

[2] We re-emphasize that the stau lifetimes at the other benchmark points would be reduced for smaller 
gravitino masses, as might be preferred by the observed light-element abundances. 
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(a) (b) 



Figure 4: T/ie tree-level Feynman diagrams for the reaction e + e — > f^fy . 

Table II: The vertices for t\tG, t^t^j and t^t^Z interactions and the corresponding tree-level 

diagrams. The parameters Cl,r depend on the stau mixing angle Of and the coefficients depend 
also on the weak mixing angle. F is the supersymmetry breaking scale. 
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IV. PRODUCTION OF STAU PAIRS IN e+e" ANNIHILATION 



The reaction e + e~ — > r+f-f proceeds via direct-channel 7 and Z exchange as shown 
in Fig. [H The tree-level vertex factors can be parametrized as given in Table fl] for the 
stau-photon-stau (t^jt^) and stau-Z-stau (r-^Zf-f) interactions. 

The cross section for the process e + e" — > r^rf with a left (L)- or right (R)-polarized e ± 
beams at a center-of-mass energy -^i is given by [27]. 
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Table III: The total cross section in pb calculated using PYTHIA with the full ISASUGRA spectrum, 
including both initial- and final- state radiation (ISR+FSR). 
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where v e = 4 sin 2 9\y — 1 and a e = — 1. P_ and P + denote the degree of polarization of the 
e~ and e + beams, respectively. The propagator terms are given by 

P s2 P s ( s ~ ™z) m 

ZZ (s - m|) 2 + m|T| ' < Z (s - m|) 2 + m\Y\ { > 

The cross sections with unpolarized electron and positron beams in the four bench- 
mark scenarios, shown in Fig. [5], exhibit a typical /3 3 threshold suppression, where (3 = 
yjl — Ami/ s. Illustrative values of the cross sections for producing f\ pairs in the bench- 
mark scenarios e, (,rj, and 6 are shown in Table HTT1 for center-of-mass energies between 0.5 
and 5 TeV. As an example, with an integrated luminosity L int = 200 fb _1 at yfs = 1 TeV, 
the benchmark points e, ( and r\ would produce 4882, 1846 and 2150 f\ pairs, respectively. 
In general, the peak cross sections are found at center-of-mass energies ^fs ~ "im^ , namely 
yfs ~ 500, 1000 and 4000 GeV for benchmarks e, r\ and (, and 9, respectively. These would 
be good energies for high-statistics studies of stau properties. The mass of the f can be 
estimated from the threshold scan and the mixing angle from the polarized cross sections. 

The unpolarized cross sections shown in Fig. [6] have non-trivial dependencies on the 
stau mixing angle 9 f . We see that the cross sections exhibit minima as functions of 6 f 
when cosOf ~ 0.6, due to the interference of 7 and Z exchange. One might hope to use 
this dependence to extract the stau mixing angle. Unfortunately, the dependencies of the 
unpolarized cross sections on 9f are minimized when cos^ ~ 0, which corresponds to 
a mainly right-handed stau and is the region expected theoretically, as discussed in the 
previous section. Moreover, there is a twofold ambiguity in the determination of 8 f , with 
the alternative value corresponding to the fi being predominantly left-handed. Considering 
the case m fl = 154.0 GeV, similar to the mass at the reference point e, and assuming 
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Figure 5: The dependences of the tree-level cross sections for f\ pair "production on the linear- collider 
center- of -mass energy, as calculated using PYTHIA and ISASUGRA. We also show the cross sections 
for producing slow staus with (3j < 0.4 in the various benchmark scenarios, also as functions 
of the center- of -mass energy of the linear e + e~ collider. The optimal energies for obtaining the 
largest numbers of these stoppable staus are different from the energies where the maxima of the 
pair-production cross sections appear. 

cos#f = 0.60, for which the cross section exhibits a minimum and hence no left/right 
ambiguity, without radiation effects we find unpolarized cross sections of 48.3 ± 0.5 pb and 
21.3 ± 0.3 pb for ^1 = 500 GeV and = 1000 GeV with L int = 2001b- 1 , respectively, 
where the errors are statistical. As seen in Fig. El these errors on the unpolarized cross- 
section measurement would not provide an accurate determination of cos 6f or of the stau 
handedness. 

One way of determining the handedness of the fi is the r decay polarization measurement, 
as mentioned previously. More detailed analysis of the measurement of tau lepton polariza- 



17, m 



The angular distribution 



tion in hadronic decay channels has been studied in 
of the pion from polarized tau lepton can be expressed as dT/Tdz = (1 + P T z)/2, where P T 
being the polarization of the r lepton. The decay fi — > tG implies P T = sin 2 9 T — cos 2 # T . 
For the points e } ( } 7] and 6 we have P T =0.902, 0.965, 0.959 and 0.999, respectively. The 
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Figure 6: The unpolarized cross sections for e + e~ — ► ?± ff in /& at y^s = 500, 1000, 3000 and 5000 
GeV as a function of the stau mixing angle Of for four of the benchmark scenarios. The expected 
statistical errors in the cross section measurement are shown for different benchmark points. We 
take the integrated luminosities Li n % = 200 fb^ 1 (for y/s = 500 GeV), 200 fb~ x (for yfs = 1000 
GeV), 400 fb- 1 (for ^ = 3000 GeV) and 1000 fir 1 (for y/s = 5000 GeV) for the ILC and CLIC 
energy options. 
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energy distributions of the pions show a maximum in which the minimum of tau energy can 
be measured. The fraction of the momenta p n ±/p T -j e t can be calculated from the measure- 
ments of pion momentum from tracker and tau-jet momentum from the calorimetric energy 
deposition. Another would be to use a polarized e ± beams. As we can see from Fig. [3, 
the polarization of the e ± beam would offer the possibility of measuring the stau mixing 
angle (and the stau mass) with much greater accuracy. Moreover, it can be used to enhance 
the signal and reduce the backgrounds. The signal is the same flavor, opposite sign taus 
with missing transverse momentum jpx- The SM backgrounds consists of r + r _ + jpx from 
the pair production of VF-bosons and Z-bosons as well as the single production Z+ /p?. 
The total cross section for W + W~ production is large compared to typical supersymmetry 
cross sections, however it proceeds via t-channel exchange of the neutrino, and hence it is 
very forward peaked. This background can be much suppressed by using a right-handed 
polarized electron beam and appropriate cuts [ij]]. Here we note that one can combine the 
measurements via threshold scan and polarized cross sections as well as the tau polarization 
to complete the determination of the SUGRA model parameters. 

In particular, as seen in Fig. [TJ, even a simple measurement of the ratio of the cross section 
with left- and right-polarized e beams would remove the ambiguity between the left- and 
right-handed stau solutions for the total unpolarized cross section measurement that was 
noted above. Indeed, measurements of the cross sections with beams could determine 
both cos Of and with interesting accuracy. We display in Fig. [8] the error bands in the 
(m^, cos 6*f ) obtained from measurements of the polarized cross sections for e + e~ — > r+ff 
at -y/i = 500, 1000, 3000 GeV. The combined experimental accuracies in the mass and mixing 
angle of the stau obtained from this Monte Carlo simulation are m fl = 154.0 ±0.7 GeV and 
A cos 6^ = 0.005, respectively. Similar analysis yield m fl = 346.0 ±2.4 GeV, A cos 6^ = 0.01 
and m fl = 327.0 ± 2.4 GeV, A cos 6^ = 0.01 using measurements at y/s = 1000 GeV and 
m ?1 = 1140 ± 13 GeV and Acos^f = 0.03 using measurements at yfs = 3000 GeV, for the 
points C, rj and 9, respectively. Here, we take the integrated luminosities L int = 200 fb _1 
for v/i = 500(or 1000) GeV and L int = 400 fb" 1 for ^1 = 3000 GeV. For a high luminosity 
option (L int = 1000 fb" 1 ) for the CLIC with yfs = 3000 GeV, the stau fi can be measured 
with less errors on the mass and mixing, namely Am f = 8 GeV and A cos 6^ = 0.02 at the 
point 9. In our calculations, we take into account the degree of polarization P_ = 90% for 
the electron (which means 90% of the electrons are left-polarized and the rest is unpolarized) 
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Figure 7: The cross section of e + e~ — > f^fT" in /& a£ = 500, 1000 and 3000 GeV as a function of 
the stau mixing angle Of, for unpolarized (U) as well as 90% left- (L) and 90% right- (R) polarized 
e~ beams, and assuming vnf x = 154, 346, 327 and 1140 GeV for the benchmark points. The red lines 
show the results for both the electron (90%) and positron (60%) polarized, while blue lines show only 
electrons (90%) are polarized. The statistical experimental errors would be similar to those shown 
in Fig. [6] 

and P + = 60% for the positron. It might be possible to reduce further the error in by 
using tracking and time-of-flight measurements, as has been considered for the LHC, but 
we see from Fig. [8] that this would be unlikely to improve significantly the determination of 
cos 6^. 



V. DETECTION OF STAU DECAYS 



Measuring stau decays using staus that stop and decay inside a linear-collider detector 



was first discussed in 



211 ] . Measurement of, e.g., the mass of the gravitino and the decay 



tau polarization require optimizing the number of staus stopped and then decaying inside 



the collider detector, which we discuss here. As in 
detector with the structure proposed in 



211 ]. we assume a future linear-collider 



281 ]. A metastable stau may stop in the hadron 
calorimeter or in the iron yoke. The amount of absorber material _R(g/cm 2 ) in the detector, 
summed along the radial direction at different longitudinal angles, is 835 to 1250 g/cm 2 
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Figure 8: The error bands for cos Of as a function of the stau mass, as obtained from measurements 
of the cross sections for e+e~ f^ff at *Js = 500,1000,3000 GeV with left (90%)- (L), right 
(90%)- (R) and unpolarized (U) e~ beams, assuming for inputs = 154,346,327 and 1140 GeV, 
as for the point e, £, tj, 0, and the value of cos Of = 0.60. Here we take the positron polarization 
(60%) opposite of electron polarization. 

for the hadronic calorimeter (HCAL) and 1810 to 2750 g/cm 2 for the magnet return yoke 
respectively. As discussed in 
m fl is given by 

logio(Pl) = r ^ un " / "" Tiy ~ ±J , (8) 



2l| . the range R of a metastable stau with given fij and mass 
[log 10 {R/m fl ) - ci] 



where c\ = 2.087 and C2 = 3.227 for steel 



Co 



3l| . The corresponding values of /?7 for staus 
stopping in different detector parts for the benchmark points are shown in Table [TV] (The 



corresponding values for some other stau masses are given in Table II of 



2lj.) As a repre- 



sentative example, we consider below the production of slow-moving staus with (3^ < 0.4, 
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Table IV: The maximal values of /?7 below which staus would be stopped in different detector parts, 
for the benchmark scenarios e,rj,£ and 9. 



m f (GeV) 


07 (HCAL) 


/?7 (Iron Yoke) 


154 


0.38-0.43 


0.48-0.55 


346 


0.30-0.34 


0.38-0.43 


327 


0.30-0.34 


0.38-0.44 


1140 


0.21-0.23 


0.26-0.30 



that might plausibly be stopped with a linear-collider detector, so that their decays could 
be observed. 

In the absence of photon radiation, the velocity of the outgoing staus in the direct pair- 
production reaction e + e~ — > f-f is simply f3 = a/1 — 4mf/s. However, the values of the 
stau velocities are altered when initial-state radiation and a realistic luminosity spectrum 
are taken into account. Fig. [9] shows the distributions of Pj expected for the fiS pair- 
produced in the various gravitino DM scenarios at the center-of-mass energies 500 GeV (the 
design for the ILC), 1 TeV and 3 TeV (the design for CLIC). The two sets of P'j distributions 
correspond to alternative luminosity spectra yielded by beam conditions designed to optimize 
the luminosity close to the nominal center-of-mass energy (left) and the total luminosity 
(right). We see that the Pj distributions exhibit substantial low-energy tails, particularly 
in the case where the total luminosity is optimized, as shown in the right panel of Fig. [9l 

We plot in Fig. [5] the cross sections the production of staus with Pj < 0.4 via direct pair- 
production e + e _ — > f^rf accompanied by initial- and final-state radiation, which opens 
up the possibility of producing slow-moving staus. We see that the optimal center-of-mass 
energies for producing slow-moving staus via direct pair-production are 330, 730, 700 and 
2500 GeV for benchmarks e,(,V and 9, respectively. We note that the center-of-mass energy 
dependencies of these cross sections for slow-moving staus are completely different from the 
total cross sections shown in Fig. [5], reflecting the increasing difficulty of radiating sufficiently 
to produce a slow-moving stau as the center-of-mass energy increases. 

Table [V] gives estimates of the numbers of staus with Pj < 0.4, as produced directly 
via e + e~ — > r-^f-f accompanied by initial- and final-state radiation at selected nominal 
center-of-mass energies. These staus might stop in a typical linear-collider detector, where 
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Figure 9: The /?7 distributions for direct stau pair-production e + e~ — ► fj fjj" a£ different center-of- 
mass energies: a) \/s = 500 GeV, b) \fs = 1000 GeV and c) \/s = 3000 GeV. The solid, dashed, 
dotted and dashed-dotted lines correspond to the gravitino DM points e, £, r\ and 9, respectively. 
These distributions were obtained from simulations of running conditions designed (left) to opti- 
mize the luminosity close to the nominal center- of-mass energy and (right) to optimize the total 
luminosity. The low-energy tails of the P"f distributions are apparent, and much larger in the right- 
hand panel. Here, we have used the CALYPSO library [3Mj to include the luminosity spectrum into 
the event generator PYTHIA. 

their decays could then be observed. We note that the numbers of stopped staus obtained 
at the optimal center-of-mass energies for benchmark points are more than an order of 
magnitude larger than the numbers that would be stopped in measurements at center-of- 
mass energies corresponding (approximately) to the peaks of the total pair-production cross 
sections, namely 500, 1000, 1000 and 3000 GeV, respectively. 

VI. INDIRECT STAU PRODUCTION 

In addition to direct stau pair-production e + e _ — > f+f-f, staus can be produced via the 
cascade decays of heavier sparticles, e.g. sleptons, — > fire, or neutralinos, Xi ~^ ^i T - 
Each heavier sparticle eventually yields one stau among its cascade decay products, which 
later decays into a gravitino. For example, for the nominal ILC center-of-mass energy 
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Table V: The numbers of staus with /3j < OA, that may be stopped in a generic linear- collider 
detector, for the different benchmark points at different center- of-mass energies and luminosities. 
In the fourth column, the energies are optimized for the direct pair-production of slow-moving staus, 
whereas in the last column all sparticle pair-production cross sections are included. 



v^GeV) 


500 


1000 


3000 


Optimal for 


Maximal including 


Wfb- 1 ) 


200 


200 


400(1000) 


pair prod'n 


other prod'n processes 


e 


34 


4 


4(10) 


1500 


1700 


c 




12 


4(10) 


254 


700 


V 




10 


4(10) 


370 


600 


e 






8(20) 


56(140) 


140(350) 



of y/s = 500 GeV, the most important contributions to stau-pair production at point e 
in fact come from the processes e + e~ — > XiXi an d e + e~ — > e^e^ L with cross sections of 
1.26 x 10 _1 pb and 7.44 x 10~ 2 / 6.76 x 10~ 2 pb, respectively. On the other hand, for the point 
6 at y/s = 3000 GeV, the main contribution comes from the pair production of right-handed 
selectrons, smuons and neutralinos with the cross sections 8.49 x 10~ 4 pb, 7.31 x 10~ 4 pb and 
9.66 x 10~ 4 pb, respectively. We display in Fig. [10] the dominant cross sections for sparticle 
pair-production in unpolarized e + e~ collisions in benchmarks e, 77, ( and 8, as functions of 
the center-of-mass energy. 

The pair-production of these heavier sparticles would be interesting in its own right, 
providing new opportunities to study sparticle spectroscopy with high precision, and to 
measure many sparticle decay modes. Some prominent decay modes of selectrons, smuons, 
neutralinos, charginos and heavier staus are shown in Table|VT] 3 : measuring these would give 
much information about sparticle masses and couplings, and hence the underlying pattern 
of supersymmetry breaking. 

Concerning the point 8, the angular distributions of the staus and right-handed selec- 
trons produced directly are given in Fig. [TT1 The staus peak around the central part of 
the distribution for their pair production while selectrons place forward direction for the 
production due to the contribution from the neutralino exchange in the t— channel. 

[3] In the case of point 9, the small — f\ mass difference implies that the only three-body decays allowed 
kinematically are £r — * f\W, and similarly for jiR. 
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Figure 10: The cross sections for pair-production of super symmetric particles in the benchmark 
scenarios e, £,7/ and 6, as functions of the e + e~ center- of -mass energy. Also shown are the yields 
of slow-moving staus with /?7 < 0.4. 

If one study only the direct production of staus it is useful to apply a cut on the angular 
distribution cos 6 < 2/3 for an easy identification of staus. In this case the direct signal 
cross section reduces 15% and the contribution from indirect production can be eliminated 
mostly. The right-handed smuons show the same angular distribution as the staus for their 
pair production. At this point, in order to be sure that stau is still an NLSP we need to 
have experimental identification of the corresponding final states. Assuming the stau is an 
NLSP we expect hadronic decay channels of tau leptons and the missing energy, otherwise 
smuons give their unique signatures in the muon spectrometer when they stopped inside the 
detector. However, here we take the stau as an NLSP with a good motivation from the BBN 
implications. We show in Fig. d0] all the contributions to the stau pairs from the indirect 
production mechanisms for the corresponding benchmark points. In addition, our interest 
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Table VI: Significant decay modes and branching ratios of selectrons, smuons, neutralinos, charginos 
and heavy staus for the chosen benchmark points e£,r) and 6. Here, the charged lepton l~ may be 
either e~ or [i~ , and the branching ratios add up to ~ 100 % when these and charge-conjugate 
modes are added. 



Sparticle 


Significant decay modes 


BR (%) for e,C,V,0 


il/h 


xii~/x^i 


100, 100, 100, 100 


~ l ~R 


ffr+r 


45,42,40,0 






55,58,60,0 






0,0,0, 100 


v T 


X\Vt 


55,21,17,64 






45, 79, 83, 36 




X\r- 


50,20,17, 56 




ffZ° 


23,37,39,16 




fffl 


27,43,44,28 




ifr* 


35,30,29,17 






7, 10,11,17 


Xi 




13,4,5,0 




v>ll~ 


18,16,16,17 




i~n 


9,13,14, 16 




V T T~ 


22,19,19,18 






8,15,15,17 


X°2 




7,2,2,0 




v T v T lv T v T 


10,10,10,9 






9,8,8,8 






5,7,7,8 






4,8,8,8 



here is in the contributions of these heavier sparticles to the total yields of slow-moving 
staus with (3^ < 0.4, which are also shown in Fig. [101 Comparing these yields with those 
produced by direct pair-production alone, shown in Fig. [5], we see that they are greatly 
enhanced. It is also apparent that the optimal center-of-mass energies are substantially 
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Figure 11: The angular distributions of the staus (solid line) and right-handed selectrons (dashed 
line) from their pair productions at the center of mass energy y/s = 3000 GeV. 

higher, reflecting the higher thresholds for pair-producing heavier sparticles. Specifically, we 
find optimal center-of-mass energies of 430, 1200, 1150 and 2800 GeV for the benchmark 
points e,(,V and 9, respectively. The last column of Table [Vl displays the maximal number 
of slow-moving and hence stoppable staus that might be obtainable at these optimal center- 
of-mass energies in the four benchmark scenarios. These numbers provide greatly improved 
prospects for observing stau decays into gravitinos. This would yield direct insight into the 
mechanism of supersymmetry breaking, via a measurement of the stau lifetime, and possibly 
of the gravitino mass, via the kinematics of stau decays. 



VII. CONCLUSION 



We have discussed in this paper how the study of a metastable stau NLSP could be 
optimized at a linear collider such as the ILC or CLIC. The stau production rate, particularly 
using a polarized e~ beam, would provide information on the handedness of the stau, as well 
as its mass. The center-of-mass energy could then be optimized for the combined direct and 
indirect production of slow-moving staus that would stop in the detector. Their decays could 
provide valuable information, not only on the stau lifetime and on the mass of the gravitino, 
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but also supplementary information on the handedness of the stau via a measurement of 
the polarization of the tau decay products. The following are possible scenarios for studying 
metastable stau scenarios at the LHC and a subsequent linear collider. 

If there is a metastable stau weighing up to ~ 300 GeV, as in benchmarks rj and (, the 
LHC has been shown able to produce and detect them, and also to analyze the mSUGRA 
models in which they appear {33I, [yj]. Specifically, it should be possible to measure the stau 
mass with an accuracy of ~ 1 %. However, most of the LHC staus would be produced with 

~ 1, and would decay outside the detector if they have lifetimes > 10~ 7 s. We recall 
that there are no significant cosmological constraints on such staus if their lifetimes < 1 s, 
and that there are only weak constraints on staus with lifetimes < 10 3 s. Therefore, collider 
information on the stau lifetime would be interesting for cosmology as well as the theory 
of supersymmetry breaking. However, the numbers of slow staus that stop inside the LHC 
detectors would be relatively modest: perhaps just a handful in benchmarks 77 and (, and 
somewhat more in benchmark e. The LHC would detect no staus in the case of benchmark 
0. 

If the stau is indeed found by the LHC, in these scenarios one would know the optimal 
energy for producing the optimal number of stopping staus via the reaction e + e~ — ► f^f^. 
As we have emphasized in this paper, there would in general be many additional staus 
coming from the decays of heavier sparticles produced by reactions such as e + e" — > e^eL, 
e + e~ — > xXi e tc The masses of these heavier sparticles could also be measured at the LHC 
in models resembling benchmarks e,T] and (, and hence the corresponding e + e~ thresholds 
could be estimated. These production mechanisms would also yield slow staus that should 
be included in the choice of the e + e~ center-of-mass energy in order to maximize the number 
of stopped staus. 

Alternatively, the stau might be too heavy to be produced at the LHC, as in benchmark 
scenario 6 which, we recall, has the added merit of being 6,7 Li-friendly. In such a case, one 
would need to measure the stau mass in e + e~ annihilation itself before going to the optimal 
energy for stopping staus. This would be depend on the masses of heavier sparticles, which 
could be determined either directly or by invoking some model. 

The signature for a stopped stau would be its decay fx — > Gr out of coincidence with 
a collision in the central detector {21]. The r decay modes (and their branching ratios) 
observable in the collider detector would be the leptonic 3-body decays r — > v T ^iv^ (17.4%) 
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and r — > v T ev e (17.8%), and hadronic decays r — > v T n (11.1%), r — > v T p — > Ur^Tr (25.4%) 
and r — > u T ?>-n — > z/ T 7r ± 7r + 7r _ + u t 'k ± 'k 'k (19.4%). In the benchmark scenarios considered 
here, the decaying rs would be very energetic, since ttiq{= 20,100,20,66 GeV) << m fl (= 
154,346,327,1140 GeV), respectively. These r decays are therefore likely to produce iso- 
lated, high-energy hadronic or electromagnetic clusters above the threshold of the HCAL 
(Eh,em > 10 GeV), hadronic showers in the yoke (Eh > 10 GeV), or energetic fis originating 
in the HCAL or yoke (E^ > 10 GeV) |21j. The main background would come from cosmic 
rays, and some rejection could be obtained by excluding decay vertices in the outermost 



detector layers and vetoing signals initiated by external muons. As pointed out in [21|, fur- 
ther excellent background discrimination would be provided by requiring the decay vertex 
to be consistent with the estimated stopping point of a f\ detected and measured previ- 
ously. The r recoil energy E T = m f (l — m^/mf) /2 would have sensitivity to the gravitino 
2l|. We note in addition that the tau decay spectra measurable in each of the above 



mass 



decay channels would be sensitive to the tau polarization and hence the handedness of the 
parent stau. 

This analysis demonstrates once again the complementarity of the LHC and a linear 
electron-positron collider of sufficiently high energy. Information from the LHC may not 
only establish the threshold for new physics, but also provide an estimate of the optimal 
energy for studying it. At least for certain studies, such as the decays of stopped staus 
considered here, this optimal energy may be significantly higher than the e + e~ production 
threshold. 
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